The Kardar-Parisi-Zhang equation in the weak noise hmit: 
Pattern formation and upper critical dimension 
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We extend the previously developed weak noise scheme, applied to the noisy Burgers equation 
in ID, to the Kardar-Parisi-Zhang equation for a growing interface in arbitrary dimensions. By 
means of the Cole-Hopf transformation we show that the growth morphology can be interpreted 
in terms of dynamically evolving textures of localized growth modes with superimposed diffusive 
modes. In the Cole-Hopf representation the growth modes are static solutions to the diffusion 
equation and the nonlinear Schrodinger equation, subsequently boosted to finite velocity by a Galilei 
transformation. We discuss the dynamics of the pattern formation and, briefly, the superimposed 
linear modes. Implementing the stochastic interpretation we discuss kinetic transitions and in 
particular the properties in the pair mode or dipole sector. We find the Hurst exponent H=(3- 
d) / (4-d) for the random walk of growth modes in the dipole sector. Finally, applying Derrick's 
theorem based on constrained minimization we show that the upper critical dimension is d=4 in the 
sense that growth modes cease to exist above this dimension. 
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I. INTRODUCTION 



The large majority of natural phenomena are char- 
acterized by being out of equilibrium. This class in- 
cludes turbulence in fluids, interface and growth prob- 
lems, chemical reactions, processes in glasses and amor- 
phous systems, biological processes, and even aspects of 
economical and sociological structures. 

As a consequence much of the focus of modern statis- 
tical physics, soft condensed matter, and biophysics has 
shifted towards such systems Q, . Drawing on the case 
of static and dynamic critical phenomena in and close to 
equilibrium, where scaling and the concept of universal- 
ity have successfully served to organize our understand- 
ing and to provide a variety of calculational tools 
a similar strategy has been advanced towards the much 
larger class of nonequilibrium phenomena with the pur- 
pose of elucidating scaling properties and more generally 
the morphology or _£attern formation in a driven nonequi- 
librium system 



field h{r,t) has the form 

dh _2, A, 

(77(r, t)rj{r', t')) = M'^(v - r')5{t - t'). 



-Vh-Vh-F + T], 



(1.1) 
(1.2) 



Here the damping coefficient or surface tension v char- 
acterizes the linear diffusion term, the parameter A con- 
trols the strength of the nonlinear growth term, is a 
constant imposed drift term, and 77 a locally correlated 
white Gaussian noise modelling the stochastic nature of 
the drive or environment; the noise correlations are char- 
acterized by the noise strength A 0, |^ 0, 0, ^| . 
In terms of the vector slope field 



Vh, 



(1.3) 



or p£ 

lii 



There is a particular interest in the scaling proper- 
ties and general morphology of nonequilibrium models 
0, 1^. Here the Kardar-Parisi-Zhang (KPZ) equation 
has played a prominent and paradigmatic role. The KPZ 
equation describes aspects of the nonequilibrium kinetic 
growth of a noise-driven interface and provides a sim- 
ple continuum model of an open driven nonlinear system 
exhibiting scaling and pattern formation [^Q- 

The KPZ equation for the time evolution of the height 
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the KPZ equation rnaps onto the Burgers equation driven 
by conserved noise 0, 0, 0, 

^ i^V^u + A(u • V)u + Vr/. (1.4) 

In the deterministic case for 77 = the Burgers equation 
has been used to study irrotational fluid motion and tur- 
bulence [ia, l2^ and also as a model for large 
scale structures in the universe plj . 

In a series of papers we have analyzed the one dimen- 
sional noisy Burgers equation for the slope field of a grow- 
ing interface. In Ref. we discussed as a prelude the 
noiseless Burgers equation 0| in terms of its non- 
linear soliton or shock wave excitations and performed 
a linear stability analysis of the superimposed diffusive 
mode spectrum. This analysis provided a heuristic pic- 
ture of the damped transient pattern formation. As a 
continuation of previous work on the continuum limit 
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of a spin representation of a solid-on-solid model for a 
growing interface [2^ , we applied in Ref . ^2^] the Martin- 
Siggia-Rose formalism [2^ in its path integral formula- 
tion 2?, "27, "2^ to the noisy Burgers equation 15, Ts"] and 
discussed in the weak noise limit the growth morphology 
and scaling properties in terms of nonlinear soliton or do- 
main wall excitations with superimposed linear diffusive 
modes. In Ref. |29j | we pursued a canonical phase space 
approach based on the weak noise saddle point approxi- 
mation to the Martin-Siggia-Rose functional or, alterna- 
tively, the Freidlin-Wentzel symplectic approach to the 
Fokker-Planck equation [sol l3ll | . This method provides 
a dynamical system theory point of view [s^ IssL Is^ l to 
weak noise stochastic processes and yields direct access to 
the probability distributions for the noisy Burgers equa- 
tion; brief accounts of the above works has appeared in 
Refs. [3^ Is^. Further work on the scaling properties 
and a numerical investigation of domain wall collisions 
has appeared in in Refs. 37, 38, 39, 40, 41] . A detailed 
summary and further developments have been given in 
an extensive paper li^ . 

In the present work we address the KPZ equation for 
a growing interface in arbitrary dimensions. Applying an 
extended form of the canonical weak noise approach in or- 
der to incorporate multiplicative noise and drawing from 
the insight gained by the analysis of the ID noisy Burgers 
equation, we identify the localized growth modes for the 
KPZ equation. The growth modes are spherically sym- 
metric and are equivalent to the domain walls or solitons 
identified in the ID case. The growth modes propagate 
and a dilute gas of modes constitute a dynamical network 
accounting for the kinetic growth of the interface. 

We also consider the issue of an upper critical dimen- 
sion for the KPZ equation. The KPZ equation lives at 
a critical point, conforms to the dynamical scaling hy- 
pothesis ^, JLQj 43J and is characterized by the scaling 
exponents z and ^|. Dynamic renormalization 

group calculations yield d = 2 as lower critical dimension 
|7, 13]. In addition to the scaling properties in the rough 
phase, characterized by a strong coupling fixed point, a 
major open problem remains the existence of an upper 
critical dimension [3 EM In the present context 

we interpret the upper critical dimension as the dimen- 
sion beyond which the growth modes cease to exist. On 
the basis of a numerical analysis and an exact argument 
based on Derrick's theorem l47| we propose that d = 4 is 
the upper critical dimension for the KPZ equation. 

The paper is organized in the following way. To bring 
the reader up to date we review in Sec.^the KPZ equa- 
tion with emphasis on the scaling properties. In Sec. IIIII 
we summarize the weak noise approach including an ex- 
tension to the case of multiplicative noise in order to 
treat the KPZ equation in the Cole-Hopf representation. 
In Sec. IIVI we address the KPZ equation and the associ- 
ated noisy Burgers and Cole-Hopf equations within the 
weak noise scheme and derive the fundamental determin- 
istic field equations governing the weak noise behavior. 
In SecElwe turn to the solutions of the field equations. 



As a prelude we review the solutions in the ID case and 
then turn to the KPZ equation in its Cole-Hopf repre- 
sentation in higher dimensions. In Sec. IVII we form a 
dynamical network of growth modes accounting for the 
growth morphology of the KPZ equation. We establish a 
field theory based on the picture of the growth modes as 
charged monopoles. Finally, we discuss briefly the super- 
imposed linear mode spectrum. In Sec. lVIll we turn to the 
stochastic interpretation and discuss kinetic transitions 
and in particular the anomalous diffusion and scaling in 
the dipole sector. In Sec. IVIIII we discuss the issue of 
the upper critical dimension and present, using Derrick's 
theorem based on constrained minimization, an algebraic 
proof of the upper critical dimension being equal to four. 
Sec. lIXl is devoted to a summary, a list of open problems, 
and a conclusion. In appendix IXI we consider the appli- 
cation of the weak noise method to Brownian motion and 
the overdamped oscillator. Aspects of the present work 
has appeared in Ref. [isl ]. 



II. THE KPZ EQUATION 

The KPZ equation (|1.1|) was proposed as a model 
for the kinetic nonequilibrium growth of an interface 
driven by noise 0, IS 0- Although the equation only 
describes limited aspects of true interface growth and, 
for example, ignores surface diffusion the equation 
has achieved an important and paradigmatic status in 
the theory of nonequilibrium processes |ll| . In many re- 
gards the KPZ equation serves as a prototype continuum 
model for nonequilibrium processes in much the same 
way as the Ginzburg-Landau functional in, for example, 
the context of critical phenomena |0i Q . 



A. General Properties 

From a structural point of view the KPZ equation Hl.l|) 
has the form of a noise-driven diffusion equation with 
a simple nonlinear term added. Whereas the diffusion 
term i/V^/i gives rise to a local flattening or relaxation 
of the interface, corresponding to a surface tension, the 
crucial nonlinear term (A/2)V/i • V/i accounts for the 
lateral growth of the interface Q. In that sense the 
KPZ equation is a genuine kinetic equation describing 
a nonequilibrium process in the sense that the drift term 
i^V^/i -I- (A/2)V/i • V/i — F cannot be derived from an 
effective free energy. The noise drives the height field 
into a stationary state whose distribution is not known 
in detail except in ID, where it is independent of A and 
given by [HEI 



Po{h) (X exp 



(2.1) 
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In the linear case for A = the KPZ equation l|l.Hl re- 
duces to the Edwards- Wilkinson equation (EW) ^50J 



V 



T = v j df^x {Vhf, (2.2) 



which in a comoving frame with velocity —F describes an 
interface in thermal equilibrium at temperature T = A 
with stationary Boltzmann distribution given by Eq. 
()2.1(l . In the EW case we also have access to the time- 
dependent distribution. Expanding the height field on a 
plane wave basis, h{r) = f d'^x exp(ikr) and introduc- 
ing the diffusive mode frequency lo\^ = i/k^ we obtain for 
the transition probability from an initial profile /i^ to a 
final profile /ik in time T |^ 



P{hk,T)^ 
A 



|/ik-<exp(-c^kT)|^ 



{27Ty 1 - exp(-2wkr) 

and, for example, the height correlation function 

A 



(M)(k,tj) = 



(2.3) 



(2.4) 



with a Lorentzian line shape controlled by the diffusive 
poles at u; — ±ivk'^] see also Ref. [5ll |. 

Averaging the KPZ equation in a state at times where 
transients have died out, we obtain 



d{h) 
dt 



{X/2)(Vh-Vh) ~ F, 



(2.5) 



showing that the nonlinear term gives rise to nonequilib- 
rium growth determined by the magnitude of (Vft. • V/i); 
note that choosing F to balance the nonlinear term, i.e., 
F — (A/2)(V/i • V/i), is equivalent to choosing a comov- 
ing frame in which d{h) /dt = 0. 

In addition to being invariant under time and space 
translations, the KPZ equation is also invariant sub- 
ject to the nonlinear Galilei transformation: 



r - 
h- 
F 



r 



u" • r, 

- (A/2)uO 



(2.6) 
(2.7) 
(2.8) 



Hence, the transformation to a moving frame with veloc- 
ity Au° is absorbed by adding a constant slope term • r 
to the height field h and shifting the constant drift term 
by (A/2)uo. Note that the invariance is associated with 
the nonlinear parameter A which enters as a structural 
constant in the Galilei group transformation; in the EW 
case this invariance is absent. 

The KPZ equation is characterized by the parameters 
J/, A, F, and A of dimension [ly] = L'^/T, [A] = L/T, 
[F] = L/T, and [A] = L'^+^/T. By transforming to a 
comoving frame ft, — > /i — _Fi we can absorb the drift 
F for a given growth morphology and from the remain- 
ing dimensionfuU parameters form the dimensionless pa- 
rameter AX'^ / 1^'^'^^ . Consequently, the weak noise limit 
A ^ is equivalent to the weak coupling limit A ^ (or 
the limit of large damping ^ oo). In Fig. ^ we have in 
2D depicted a growing interface. 
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FIG. 1: We depict a growing interface in 2D (arbitrary units). 



B. Burgers equation and Cole-Hopf equation 

The paradigmatic importance of the KPZ equation 
stems from the fact that it in addition to describing 
nonequilibrium surface growth also is associated with the 
theory of turbulence via its equivalence to the noisy Burg- 
ers equation H1.4|l . Moreover, by means of the nonlinear 
Colc-Hopf tranformation |2|, |^ |^ 



w = exp 



2v 



(2.9) 



the KPZ equation takes the form of a linear diffusion 
equation driven by multiplicative noise, denoted here the 
Colc-Hopf equation, 

— ^uV^w wFH wrj. (2.10) 

dt 2v 2v 

In the absence of noise Eq. H2.10|l becomes the linear 
diffusion equation and is readily analyzed, thus allowing 
a discussion of the KPZ and Burgers equations in the 
deterministic case 0, • In the noisy case a path in- 
tegral representation of the solution of Eq. (|2.1(J|) maps 
the Cole-Hopf equation and thus the KPZ equation onto 
a model of a directed polymer with line tension 1/4^ in 
a quenched random potential (\/2v)r](y,t^. The disor- 
dered directed polymer model is a toy model within the 
spin glass literature and has been analyzed by means of 
the replica method and Bethe ansatz techniques 11, "531. 

The Galilean invariance of the KPZ equation is for the 
Burgers equation H1.4|l and the Cole-Hopf equation (|2.1U|) 
supplemented by the transformations 



u -> u-f u°, 

w ^ w exp [(A/2i/) u° • r] , 



(2.11) 
(2.12) 



4 



i.e., the slope field is shifted by u" and the Cole-Hopf 
field by the multiplicative factor exp[(A/2z/)u'' • r]. 



C. Scaling properties 

The KPZ equation conforms to the dynamical scaling 
hypothesis with long time-large distance height correla- 
tions HEi 



{hh){r,t) :^r^'^^t/r'-) 



(2.13) 



characterized by the roughness exponent dynami- 
cal exponent z, and scaling function $. Consequently, 
most theoretical efforts have addressed the scaling issues. 
Based on i) perturbative dynamic renormalization group 
(DRG) calculations in combination with the scaling law 



C + z = 2, 



(2.14) 



following from Galilean invariance, and the known sta- 
tionary height distribution in ID given by Eq. H2.ll) 
0, El Elli ii) the mapping of the KPZ equation onto 
directed polymers in a quenched environment and ensu- 
ing replica calculations [ll|, |5J| , iii) mode coupling calcu- 
lations [53, m, m , iv) operator expansion methods 
[6lj | , and v) numerical calculations [gJ 0, 113 , the 
following picture has emerged. 

In ID the interface is rough and characterized by a per- 
turbatively inaccessible strong coupling fixed point with 
scaling exponents z — 3/2 and ( — 1/2, following from 
the stationary distribution in combination with the scal- 
ing law. Above the lower critical dimension d = 2 a 
DRG calculation in e = d — 2 predicts a kinetic tran- 
sition line between a smooth phase characterized by a 
weak coupling fixed point at A = with exponents z = 2 
and C = (2 — d)/2 and a rough phase characterized by a 
poorly understood strong couplin g fix ed point. On the 
transition line z — 2 and C = |44l66l| . Based on numer- 
ics the following expressions for z have been proposed: 
Zkk = 2{d + 2)/{d + 3) 67] and zwk = (2d+ l)/(d+ 1) 
[68l| . Both Zkk and zwk agree with z = 3/2 in d = 1; 
for d 00 we have z 2, corresponding to a smooth 
phase at an infinite upper critical dimension. An oper- 
ator expansion method l6lj yields Zl = id/ {2d + 1) for 
d > 1 and z ^ 2 for d — > cx3. The scaling properties of 
the KPZ equation is summarized in Fig. |21 in a plot of 
the renormalized coupling strength g = X^A/v^ versus 
spatial dimension. 

Whereas d — 2 constitutes the lower critical dimen- 
sion permitting a loop expansion in powers of e, the is- 
sue of an upper critical dimension has been much de- 
bated. Mode coupling approaches yield an upper critical 
dimension d = 4 with a possible glassy behavior above 
4 m m m IToI. loop expansion to all orders in e 
supported by an exact evaluation of the beta-function in 
a Callen-Symanzik renormalization group scheme asso- 
ciates the uppe r critical dimension with a mathematical 
singularity 0,E3- 




STRONG COUPLING 
FIXED POINT 



(d>1) 



KINETIC PHASE 
TRANSITION 



WEAK 
(SMOOTH) 



^ = 
I 



FIG. 2: We summarize the scaling properties of the KPZ 
equation in a plot of the renormalized coupling constant (the 
fixed point value) AA'^/z/^ versus the dimension d of the sys- 
tem. Below d = 2 the scaling properties are determined by 
the strong coupling fixed. Above d = 2 the system exhibits 
a kinetic phase transition from a weak coupling phase to a 
strong coupling phase. We have also indicated three conjec- 
tures for the exponent z: The operator expansion conjecture 
zl by Lassig |6J| and the conjectures based on numerics, zwk 
by Wolf and Kertesz (g^ and zkk by Kim and Kosterlitz |6^ . 



Two outstanding issues thus remain as regard the scal- 
ing properties of the KPZ equation: The upper critical 
dimension and the properties of the strong coupling fixed 
point above d = 1, i.e., the scaling exponents and scaling 
function. There is, moreover, the more general question 
of the deeper mechanism behind the stochastic growth 
morphology. In that respect the KPZ equation shares its 
strong coupling features with another notable problem in 
theoretical physics: Turbulence. 

Since the DRG in its perturbative form as an expan- 
sion about the linear theory fails to yield insight into 
the strong coupling features there clearly is a need for 
alternative systematic methods. Here we should like 
to emphasize that both the mode coupling approaches 
m EE m m m and operator expansion meth- 
ods [61| do not quahfy as being systematic. The mode 
coupling approach is based on a truncation procedure 
ignoring vertex corrections and thus violating Galilean 
invariance which is essential in delimiting the KPZ uni- 
versality class. The operator product expansion imposes 
an ad hoc mathematically motivated operator structure. 



III. 



THE WEAK NOISE METHOD 



The weak noise method is based on an asymptotic 
weak noise approximation to a general Langevin equa- 
tion driven by Gaussian white noise. The method dates 
back to Onsager 0, [t^ and has since reappeared as the 
Freidlin-Wentzel theory of large deviations |33,|33| and as 
the weak noise saddle point approximation to the func- 
tional Martin-Siggia-Rose scheme jH [H HHH]. The 
method known also as the eikonal approximation has. 
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moreover, been used in the context of thermally acti- 
vated escape [H [zl III [zl IH IzHzl III The 

weak noise or canonical phase space method has been 
discussed in Ref. "2^. In the present context we need 
a generalization of the method in order to accommodate 
multiplicative noise and therefore review and extend the 
method below. 



A. General properties 



The point of departure is the generic Langevin equa- 



with Hamiltonian 



tions for a set of stochastic variables w„ 
driven by white Gaussian noise j8l| | 



dwp _ _1 p 
dt ~ 2 " 



(3.2) 



where Fn{wq) is the drift, Gnm{wq) is accounting for mul- 
tiplicative noise, V„ = d /dwn, and A is the explicit noise 
strength; sums are performed over repeated indices. In 
the case of multiplicative noise with G depending on w 
this is the Stratonovich formulation. In the Ito formula- 
tion the compensating drift term GVG is absent corre- 
sponding to the application of the Ito differentiation rules 
[s^ l83l|. Note that in the case of i) additive noise, i.e., 
G independent of w or ii) to leading order in A the Ito 
and Stranovich formulations are equivalent. The generic 
Langevin equation (|3.1|) driven by white noise encom- 
passes with appropriate choices of F and G all white 
noise driven continuous Markov processes. 

On the deterministic level the corresponding Fokker- 
Planck equation for the probability distribution P(w„, t) 
has the form ^ll 



dP 1 



V„[F„ + AV™is:„„]P, 



dt 2 

where the symmetrical noise matrix is given by 

Kpm{Wq) — Gpn{Wq)Gmn{Wq). 



(3.3) 



(3.4) 



In complete analogy with the WKB approach in quan- 
tum mechanics js^l, ^ oc exp[iS'/fi,], relating the wave 
function to the classical action S evaluated on the basis 
of the classical Hamilton equations of motion, it is useful 
to capture weak noise effects and relate the stochastic 
problem to a scheme based on classical equations of mo- 
tion by means of a weak noise WKB approximation to 
the Fokker-Planck equation H3.3|l . Thus introducing the 
Wentzel-Kramers-Brillouin (WKB) ansatz 



P(w„,t) cx 



S{Wn,t) 

A 



(3.5) 



we obtain to leading order in the noise strength A the 
Hamilton- Jacobi equation 0, [s^ 



_9 
di 



+ H = , p„ = VnS, 



(3.6) 



1 1 



(3.7) 



The Hamiltonian equations of motion follow from 
dwn/dt = dH/dpn and dpn/dt = —dH/dwn, 



^-A„ F +K V 



(3.8) 



^Pm^ nFm ^PmPq^ n^rnq: (^■'^) 

determining classical orbits on the energy surfaces given 
by 7J in a classical phase space {{wn}-, {Pn})- Finally, the 
action S is given by 



S{Wn,T) 



dt p„ 



dWn 

IT 



HT, (3.10) 



yielding according to Eq. H3.5|l the transition probability. 
By means of the equation of motion (|3.8I) the action can 
be reduced to the form 

S{W„,T)^-I dt KnrnPnPm- (3.11) 

The weak noise scheme bears the same relationship 
to stochastic fluctuations as the WKB approximation in 
quantum mechanics, associating the phase of the wave 
function with the action of the corresponding classical 
orbit [s^l . In addition to providing a classical orbit pic- 
ture of stochastic fluctuations and thus allowing the use 
of dynamical system theory ^2, i54, M., EL M\ , the 
method also yields the Arrhenius factor P oc exp(— S'/A) 
for a kinetic transition from wjj to Wn during the transi- 
tion time T . Here the action S serves as the weight in the 
same manner as the energy E in the Boltzmann factor 
P oc exp(— /3i?), (3 = l/fcP, for equilibrium processes. 

In the weak noise scheme the stochastic Langevin equa- 
tion l|3.1|l is replaced by the deterministic equation of 
motion H3.8|) for w„ which together with the equation of 
motion (|3.9|l for the canonically conjugate noise variable 
Pn determine orbits lying on the constant energy mani- 
folds H = const, in a canonical phase space spanned by 
Wn and Pn- Determining a specific orbit from to w„ in 
time T by solving the equations of motion (|3.8|) and 13.9|l 
with slaved variable, evaluating the action H3.10|l . 

then yields the (unnormalized) transition probability in 
Eq. (TTKll . 

A stationary distribution of the Fokker-Planck equa- 
tion (|3.3|l is given by 



Po{Wn 



lim P{wn,T), 



(3.12) 



satisfying (1/2)V„[F„ + AV„X„„,]P = 0. Within the 
dynamical system theory framework to leading asymp- 
totic order in A this implies that S{wn,T) — s- S'o(iy„) for 
fixed final configuration w„, i.e., Po oc exp[— S'o/A]. To 
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FIG. 3: The {{w„}, {pn}) phase space with a finite-time orbit 
firom {wJi} to {w„} in time T on a _ff 7^ energy manifold, 
yielding the transition probability, and an infinite-time sta- 
tionary orbit to {wn} passing through the saddle point SP 
and lying on the H — energy manifold. 



is a weak noise approximation the scheme, like WKB in 
quantum mechanics, is nonperturbative and is thus ca- 
pable of accounting at least qualitatively for strong noise 
effects. This point of view was stressed by Coleman in 
the context of quantum field theories 0| . 

This completes our general discussion of the weak noise 
approach. As an illustration we consider in appendix IXI 
the application of the method to two systems with a sin- 
gle degree of freedom: (i) random walk and (ii) the over- 
damped oscillator. See also an application to a nonlin- 
ear finite-time-singularity model in Refs. [oil l93 | and to 
an extended system, the noise-driven Ginzburg-Landau 
model, in Refs f93L l94lj . 

IV. THE KPZ EQUATION FOR WEAK NOISE 

Here we apply the weak noise method to the KPZ equa- 
tion, the corresponding noisy Burgers equation and Cole- 
Hopf equation. 



attain the stationary limit an orbit from an initial config- 
uration wjj to a final configuration Wn traversed in time 
T — !• 00 must (i) pass through a saddle point and (ii) 
in the limit lie on a zero-energy manifold H = 0. The 
zero-energy manifold is in general composed of two sub- 
manifolds intersecting at the saddle point. The transient 
or noiseless submanifold p = 0, yielding H — and con- 
sistent with the equations of motion, corresponds to the 
transient motion determined by the noiseless equation of 
motion dwn/dt = — (l/2)p„i^„. The stationary or noisy 
submanifold determined by the orthogonality condition 
Pn[KnmPm — Fn] = 0, corrcspouds to the stationary mo- 
tion. The loss of memory and Markovian behavior result 
from the infinite waiting time at the saddle point; for 
details see e.g. Ref. poj . 

We, moreover, note that the weak noise scheme has a 
symplectic structure with Poisson bracket algebra [85ll89| 

{Wn,P„i} ^ Snm, {^n , Wm} = {Pn , Pm } = 0, (3.13) 



and Hamiltonian equations of motion 

dWn r . dpn . , 

— = {w„,i/}, — ={p„,i/}. 



(3.14) 



In Fig. 13 we have depicted the generic phase space struc- 
ture in the case where the system possesses a stationary 
state. 

The purpose of the weak noise scheme is twofold. On 
the one hand, it provides an alternative way of discussing 
stochastic phenomena in terms of an equivalent deter- 
ministic scheme based on canonical equations of motion, 
orbits in phase space, and the ensuing dynamical system 
theory concepts. The scheme, on the other hand, also 
provides a calculational tool in determining the transi- 
tion probabilities and ensuing correlations. It is also im- 
portant to keep in mind that although the starting point 



A. The KPZ case 

Adapting the KPZ equation (|l.l(l to the weak 
noise scheme we make the assignment: w„(t) 
h{r,t),p,,{t) ^ Pir, t),Knm ^ '^'^(r - r'),F„ ^ 
-2 [i^V^h + {X/2)Vh -Vh- F]. Here the index n be- 
comes the spatial coordinate r. We thus obtain the KPZ 
Hamiltonian density 



n 



KPZ 



i^VVi + -Vh -Vh- F+-P 



and the canonical field equations 



— = „V^h + -Vh -Vh-F + p, 
at 2 



dp 



-vV'^p + \Vp ■ Vh + ApV^/i. 



(4.1) 

(4.2) 
(4.3) 



The reduced action and transition probability are given 
by 



^KPZ/ 



1 



{h,T) = - d^xdt pir,t)\ (4.4) 



P^^^(/i,r) ocexp 



5^^^(/i,r) 



(4.5) 



The height field h and noise field p are canonically conju- 
gate variables in the canonical phase space ({h{r}, {p{r}), 



{h{r),p{r')}^d''{r~v'). 



(4.6) 



We, moreover, have the generators of time translation 
and space translation, i.e., energy and momentum. 



dh 



(4.7) 



jjKPZ ^ J ^ ^j^^ jjKPZ} = Vh. (4., 
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The Galilean invariance of the KPZ equation Hl.l|l im- 
plies that the noise field p is invariant; this also follows 
from the invariance of the action in Eq. 14.4|l . 

This completes the formal application of the weak noise 
scheme to the KPZ equations. The resulting classical 
field theory must then be addressed in order to eventually 
evaluate transition probabilities and correlations. 



B. The Cole-Hopf case 

The Cole-Hopf equation (|2.1()() is obtained by ap- 
plying the nonlinear Cole-Hopf transformation H2.9|) to 
the KPZ equation Hl.l|l . The Cole-Hopf equation is 
driven by multiplicative noise and most work has been 
based on the mapping to directed polymers in a ran- 
dom medium 11]. In the present context it turns 
out that a weak noise representation provides a par- 
ticular symmetric formulation. Hence, making the as- 
signment, Wnit) w{r,t),Pn{t) p{r,t),Knm ^ 

{X/2iyyw'^S'^{r - r'), F„ ^ -2 [i^V^w - {X/2i^)wF] , we 
obtain the Cole-Hopf Hamiltonian density 



CH 



(4.9) 



where we have introduced two characteristic inverse 
length scales 

The field equations are given by 
dw 
'dt 

dt 

and the reduced action and distribution by 



j'[V^w — k'^w] + k.QW'^p, 
i^[V^p — k'^p] — k^p^w, 



(4.10) 

(4.11) 
(4.12) 



1 



w,T 



{w,T) (X exp 



dtd'^x (w(r,i)p(r,i))2,(4.13) 
"'"'•"•^'1 (4.14) 



The Cole-Hopf and noise fields w and p, spanning the 
canonical phase space ({w(r)}, {p(r)}), satisfy the Pois- 
son bracket 

Mr),p(r')} = <5'^(r-r'), (4.15) 
For the total energy H'^^ and momentum 11'^^ we have 

i/™ = J d'^x , - ^, (4.16) 

U'^^ ^ J d'^x wVp , {w,n^"} = Vw. (4.17) 

Finally, since the action H4.13|l is invariant under the 
Galilean transformations (|2.6|l and (|2.12() the noise field 
p must transform according to 

p pexp [-(A/2iy)u° • r] . (4.18) 



C. The Burgers case 

The noisy Burgers equation (|1.4|) for the noise-driven 
slope of an interface will also be needed in our weak noise 
analysis. In this case we choose the assignment: Wn{t) 
u(r,i),p„(t) ^ p(r,t),i^„™ ^ V25'^(r - r'),F„ ^ 
— 2[i/V^u -I- A(u • V)u]. Hence, the Hamiltonian density 
is 



i/V^u + A(u • V)u - V(V • p) 



, (4.19) 



and the ensuing field equations 
d 



Au- Vj u = z/V^u- V(V • u), (4.20) 
Au • V ) p = 



at 

d_ 

di 

z.V2p + A(p(V-u)-(p.V)u), 



(4.21) 



where we have used the property that u is a longitudinal 
vector field. Since the operator {d/dt~ Au • V) is invari- 
ant under the Galilei transformations (|2.6(l and (|2.1HI 
the equations of motion H4.2U|I and (|4.21|) are manifestly 
Galilean invariant with an invariant noise field p. 
For the reduced action and distribution we have 



5^(u,T) = - 



P^(u,r) = exp 



dtd'^x p(r, tf, 

s^{u,Ty 



(4.22) 
(4.23) 



For the fields u and p, spanning the canonical phase 
space ({u(r)}, {p(r)}), we obtain the Poisson bracket 



{M„(r),p„j(r')} = Sn„iS'^{r - r'). 



(4.24) 



Moreover, the total energy and momentum 11^ are 
given by 

H^ = J d^x n'' , {u, i/B} = (4.25) 

= y" d'-x UnVp„ , {u„,nC"} = Vm„. (4.26) 

We note that the action H4.22|l also implies that the noise 
field p is invariant under the Galilean transformation. 



D. Phase space behavior and canonical 
transformations 



In the KPZ, Cole-Hopf, and Burgers cases the orbits or 
trajectories from an initial configuration (/i',w',u') to a 
final configuration (/i, w, u) traversed in time T, yielding 
the actions and transition probabilities, live in the corre- 
sponding phase space spanned by the canonically conju- 
gate fields, i.e., the height field h, the diffusive Cole-Hopf 
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field ui, the Burgers slope field u and their associated 
noise fields. The canonical field theories are conserved 
and the orbits are confined to constant energy manifolds, 
TJKPZ ^ const., iJ^" = const., and = const. 

The structure of phase space determines the nature of 
the underlying stochastic model. Here the zero-energy 
manifolds, H^^^ = H'^^ = = 0, play a central role 
in determining the stationary stochastic state. First we 
notice that in all three cases a vanishing noise field p — 0, 
p = 0, and p = is consistent with the field equations 
(Oil . (lOll . (EH), Km . (11201), and fOTTl . and yield 
fjKPZ ^ fjCR = f{B ^ Q ^Yie transient zero-noise 
submanifold we thus obtain the deterministic damped 
evolution equations 



dh „ A 

— = uV^h + -Vh -Vh^F, 
at 2 

dw 



dt 
du 

'dt 



i^V^u + A(u • V)u, 



(4.27) 
(4.28) 
(4.29) 



describing the transient relaxation of the height, Cole- 
Hopf, and Burgers fields subject to transient pattern for- 
mation. On the transient zero-energy submanifold w van- 
ishes in the long wave length limit as w cx exp[— i^fc^i], 
corresponding to h = {2i'/X)\ogw ^ —Ft, i.e., setting 
the drift = by the transformation h ^ h — Ft, the 
height field likewise vanishes as does the slope field u. 

The linear Cole-Hopf equation (|4.28() is exhausted by 
diffusive modes, i.e., 



w{r,t) 



-vk-'t 



df^p 
(2^ 



A(p)e 



-up t ipr 



(4.30) 



where ^(p) is chosen according to the imposed initial 
and boundary conditions. Together with the relations 
h = (l/fco)logw and u = {l/ko)'Vw/w this constitutes 
a complete solution of the deterministic transient case. 

In ID the transient pattern formation in the slope 
field is composed of propagating right hand domain walls 
connected by ramps with super imposed linear diffusive 
modes. For the height profile this corresponds to a pat- 
tern of upward pointing cusps (domain walls) connected 
by parabolic segments (ramps) For further analysis 
in higher dimensions see e.g. Ref. [T5l |. 

In the presence of noise, corresponding to the coupling 
to the noise fields (p, p, or p), the orbits in phase space 
from an initial configuration at time i = to a final 
configuration at time T with the noise field as a slaved 
variable, veers away from the transient manifold, pass 
close to a saddle point, and asymptotically approach an- 
other zero-energy submanifold intersecting the transient 
submanifold at the saddle point. This behavior follows 
from the general discussion in Sec. lIIII and is shown in the 
case of the overdamped oscillator discussed in appendix 

El 

The noisy or stationary zero-energy submanifold de- 
termines the stationary stochastic state. In the limit 



T ^ oo the orbit from fixed initial to final configurations 
migrate to the zero-energy manifold, moves along the 
transient submanifold, passes through the saddle point 
experiencing a long (infinite) waiting time, and eventu- 
ally moves along the stationary submanifold to the final 
configuration. The infinite waiting time at the saddle 
point thus ensures Markovian behavior, i.e., loss of mem- 
ory, in that the stationary distribution only depends on 
the final configuration in the limit T — > oo. 

The stationary submanifold is determined by the or- 
thogonality condition, for example, in the KPZ case. 



H 



KPZ 



d'^xp 



-Vh-Vh-F 



0,(4.31) 



i.e., the manifold vV'^h+ {\/2)Vh -Vh- F + {l/2)p or- 
thogonal to p. In the case of the overdamped oscillator, 
discussed in appendix ^ for one degree of freedom de- 
fined by Eq. (|A8I) . H — {l/2)p{p— 2jx), the zero-energy 
submanifolds are p = and p = 2^x. The action with 
H = then immediately yields S = J dt pdx/dt = 72:^ 
and hence the stationary distribution in Eq. (jA13p . On 
the other hand, in the field theoretical case characterized 
by an infinite numbers of degrees of freedom it is in gen- 
eral a difficult task to determine the manifold p = F{h) 
satisfying the orthogonality condition H4.31|l . 

As discussed in Sec. ^ the situation is special in ID. 
Here the stationary Fokker-Planck equation in the KPZ 
case admits an explicit solution given by Eq. 112. 1|) 
[4^. Within the weak noise scheme the existence of 
this fluctuation-dissipation theorem is tantamount to 
the explicit determination of the stationary submani- 
fold. This is most easily done in the Burgers case with 
Hamiltonian density Ti.^ = p{vV'^u + \uVu - (1/2) V^p). 
Setting V^p = 2v\/'^u we have TL^ = [2vX/2,)y {u^), 
i.e., a total derivative yielding = with vanishing 
slope boundary conditions. For the action we obtain 
S = j dtdx2vudu/dt — v J dxu^ — v J dxiVh)^ yielding 
the stationary distribution in Eq. H2.1|) . see also Refs. 
J,9, 42]. 

The stochastic KPZ, Cole-Hopf, and Burgers equa- 
tions Hl.l|) . H2.10|l . and H1.4|l all constitute equivalent de- 
scriptions of a growing interface. Within the weak noise 
approach the canonical structure implies that the three 
equivalent descriptions are related by canonical transfor- 
mations. By inspection we find that the Cole-Hopf and 
KPZ formulations are connected by the canonical trans- 
formations il Ei 



w = exp[fco/i] , p = /cq ^pexp[— /cq/i], (4.32) 



together with the inverse transformations 
h ^ log w , p ^ kopw; 



(4.33) 



the generating function is Gi{p,h) — pexp[A:o/i], i.e., 
dGi — wdp -\- pdh, implying w — dGi/dp and p = 
dGi/dh. Likewise, the KPZ and Burgers formulations 
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are related by means of the transformations 

u = V/i, p=-V-p, (4.34) 

with generating function G2{h,p) = p • V/i, i.e., dG2 = 
pdh + Vft. • dp, implying p = dG2/dh, and u = dG2/dp. 

From the field equations in the Burgers case, Eqs. 
(|4.2Q(I and (|4.21() . it follows that only the longitudinal 
component of p couples to the slope field u. Since Eq. 
H4.21(l is linear in p we can without loss of generality 
just keep the longitudinal component of p, i.e., p = V0, 
where (/) is a scalar potential. We thus obtain the canon- 
ical transformation u = V/i and p = — and the 
field equation (|4.2()|l takes the form [d/dt — Au • V)u = 
i/V^u + V^(/). 



V. FIELD EQUATIONS 




FIG. 4: We show the phase space structure in the ID Burgers 
case. The transient E = G, p = Q submanifold and the sta- 
tionary E = Q, p = 2vu submanifold intersect at the saddle 
point SP. We depict a finite time orbit from u to u" in time 
T. 



The canonically conjugate field equations of motion 
in the KPZ, Cole-Hopf, or Burgers formulations form 
the central starting point for an analysis of the pattern 
formation and scaling properties of the KPZ equation. 
As discussed above the three formulations are related 
by canonical transformations and represent the same 
stochastic problem in the weak noise limit. 



A. General properties 

The field equations, determining orbits in the corre- 
sponding multi-dimensional phase space, have the form 
of coupled nonlinear hydrodynamical equations. Com- 
mon to all three formulations is that the field equations 
for the noise fields or p have negative diffusion co- 
efficients. This corresponds to a Fourier mode of the 
noise field growing exponentially in time, rende ring a nu- 
merical integration forward in time unfeasible [33 • The 
growth of the noise field is consistent with the property 
that the noise drives the system into a stationary state at 
long times corresponding to orbits in phase space leaving 
the transient submanifold, see e.g. the discussion of the 
overdamped oscillator in appendix 1X1 

Leaving aside the possibility that the coupled field 
equations in the KPZ, Cole-Hopf, or Burgers formula- 
tions are exactly integrable in the sense that a Lax pair 
and an inverse scatte ring tr ansformation can be identi- 
fied, see e.g. Ref. [l9l l88ll95l | . we note that the equations 
of motion all are invariant subject to a Galilei transfor- 
mation combined with a rescaling or shift of the fields. 
This property suggest the possibility of constructing lo- 
calized propagating solitons or elementary excitations by 
first finding static localized solutions which subsequently 
are boosted by a Galilean transformation. 



B. The one dimensional Burgers case 

In the ID case summarized in detail in Ref. the 
scalar slope field u — dh/dx turns out to be the conve- 
nient variable. Since this case has served as a theoretical 
laboratory we review it here. In the ID Burgers case the 
Hamiltonian takes the form 

= J dx p[iy\/^u + Xu\7u- {l/2)\/^p], (5.1) 
yielding the equations of motion 

- AuV^ u = i/V^w - V^p, (5.2) 

note that the last term in Eq. H4.21|l vanishes in ID. In 
Fig. ^ we have depicted the Burgers phase space struc- 
ture. 



1. Domain wall solutions and pattern formation in ID 

In addition to superimposed extended phase-shifted 
diffusive modes with dispersion uj = ivp^, the equations 
H5.2(l and l|5.3|l support two localized distinctive soliton 
or domain wall modes, in the static case of the kink-like 
form, 

k 

u(x) = ±— tanhfca;, (5-4) 
ko 

with inverse scales k and fco given by Eq. H4.10|l . Boost- 
ing a pair of well-separated matched right and left hand 
domain walls by means of the Galilean transformation 
(|2.ti|) and (|2.11|) to the velocity v = 2vk, corresponding 
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FIG. 5: We show the fundamental static right hand and left 
hand domain wall or soliton solutions for the ID Burgers 
equation. The right hand domain wall is a solution in the 
deterministic case; the left hand domain wall is induced by 
the noise. 





L U 











u 




FIG. 6: We show the propagating two-domain wall quasi par- 
ticle or excitation for the ID Burgers equation. 



FIG. 7: We show a multi-domain wall representation of a 
growing interface in the slope field and the associated height 
field. 



to the slope shift vP — fc/fco, we obtain a propagating 
elementary excitation or quasi particle with vanishing u 
at infinity. The moving domain wall pair is equivalent 
to a moving step in the height profile, corresponding to 
adding a layer to the growing interface at each passage 
of the quasi particle. In Fig. we have shown the fun- 
damental static right hand and left hand domain wall 
solutions. In Fig. we have depicted the two-domain 
wall quasi particle. 

A general growth morphology is obtained by match- 
ing a dilute gas of propagating domain walls in terms of 
the slope field; for the height field this morphology cor- 
responds to a growing interface. Superimposed on the 
domain wall pattern are extended diffusive modes. In 
the limit of vanishing nonlinearity for A = the domain 
wall gas vanishes, the growth ceases, and the diffusive 
excitations exhaust the mode spectrum. 

In terms of the height field h the static domain wall 
solutions H5.4|l corresponds to the profiles 



h{x) — ±-5- In I coshfcxl, 
fco 



(5.5) 



i.e., concave and convex cusps. The corresponding Cole- 
Hopf field w is given by 

w{x) — cosh^^ kx, (5.6) 

corresponding to a localized bound state of width 1/fc 
falling off like w{x) oc exp(— fc|x|) and a concave profile 
increasing like w(x) cx exp(A:|a;|). In Fig. we have de- 
picted a multi-domain wall representation of a growing 
interface and the associated height profile. 



2. Scaling properties in ID 

The scaling properties follow from the dynamics of the 
domain walls. The right hand domain wall corresponds 
to vanishing noise field p — and is the well-known 
viscosity-smoothed shock wave of the noiseless determin- 
istic Burgers equation _16j 22j. The left hand domain 
wall, on the other hand, is associated with the noise 
field p ~ 2vu and carries action S = {l/&)v\u^T ^ en- 
ergy E = — (1/6)i/Am'^, and momentum 11 = 0, where 
u = 2fc/fco is the domain wall amplitude. For the 
quasi particle composed of a right hand and left hand 
domain wall moving with velocity v — 2vk we obtain 
S ^ {l/&)v\u^T, E = -{2/?,)u\u^, and H = vu^ . Elim- 
inating the amplitude the quasi particle is characterized 
by the gapless dispersion law 



E ■ 



4 A 



31/1/2 



in|3/2 



(5.7) 



with power law exponent 3/2. Note that the diffusive 
modes have the dispersion law to cx vk^ , corresponding 
to the power law exponent 2. 

The scaling exponents z = 3/2 and C = 1/2 follow from 

i) a spectral representation of the slope correlations and 

ii) the structure of the zero energy manifolds, i) Drawing 
on the analogy with a quantum system with Planck con- 
stant A we invoke heuristically a spectral representation 
for the slope correlations 



{uu){x,t) = / <mF{Ti)e 



■-t+illx 



(5.8) 
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where F(n) is an appropriate form factor. Within the 
single quasi particle sector, inserting the dispersion law 
H5.7|l . it readily follows that t scales with x^^^ and we 
infer the dynamic scaling exponent z — 3/2. ii) From 
the Hamiltonian (|5.1|l we infer the zero-energy mani- 
folds p = and p ~ 2vu, consistent with the equa- 
tions of motion (|5.2|) and H5.3|) . At long times the or- 
bits in the phase space approach the zero-energy 
manifold p = 2vu and we obtain for the action S ^ 
J dxdtpdu/dt — 2v j dxdtudu/dt = v J dxu^,i.e., Eq. 
(|2.1|) . and the independent stationary fluctuations of u 
are given by a Gaussian distribution. This in turn im- 
plies that the height field h = J dx'u performs random 
walk yielding according to the scaling form H2.13|l the 
roughness exponent ^ — 1/2. Note that the scaling law 
z + = 2 is automatically obeyed since the weak noise 
formulation is consistently Galilean invariant. 

The notion of universality classes is here associated 
with the dominant gapless quasi particle dispersion law. 
The scaling properties follow from the low frequency 
(long time) - small wavenumber (large distances) limit. 
For A = there is no growth, the mode spectrum is 
exhausted by extended diffusive modes with gapless dis- 
persion uj = vk'^ , yielding according to the spectral form 
the dynamic exponent z = 2. The roughness exponent 
C = 1/2 and the scaling law is not operational. This con- 
stitutes the Edwards- Wilkinson universality class. For 
A 7^ localized domain wall growth modes nucleate 
out of the diffusive mode continuum with dispersion 
E (X (A/z/i/^)n3/2; this is the KPZ universahty class. 

Summarizing, in ID the growing interface problem can 
be analyzed in some detail and a consistent interpretation 
in the WKB sense can be advanced within the weak noise 
formulation. The approach yields: i) a many body de- 
scription of a growing interface in terms of a ID matched 
network of moving domain walls with superimposed diffu- 
sive modes, ii) scaling properties and scaling exponents 
follow from the dispersion law of the dominant gapless 
domain wall excitations and the structure of the zero- 
energy manifold in phase space, iii) universality classes 
are associated with the class of gapless excitations gov- 
erning the dynamics of the interface. In Fig. |Slwe have 
in a log-log plot depicted the domain wall and diffusive 
mode dispersion laws. 



C. The Cole-Hopf case 

In higher dimension we must address the field equa- 
tions of motion in either the KPZ formulation, Eqs. 
(14. 2() and (|4.3|l . the Cole-Hopf formulation, Eqs. (|4.11() 
and (|4.12l) . or the Burgers formulation, Eqs. (|4.2U|) and 
(|4.21l) . Based on the working hypothesis that the growth 
morphology is associated with a network of growth modes 
and drawing from the insight gained in ID, the program 
is again to search for localized solutions of the equations 
of motion. Whereas both the KPZ and Burgers formula- 
tions do not easily yield to analysis, the symmetric Cole- 



log E, log CO 




wall mode 



* log n, log k 



FIG. 8: We show in a log- log plot the domain wall dispersion 
law E oc Yl^^^ and the diffusive mode dispersion law lo (x . 

Hopf formulation turns out to be the convenient starting 
point 

1. Static localized modes 

In the static limit the Cole-Hopf equations of motion 
(|4.11(l and H4.12|l assume the symmetrical form 

j/V^w = vk^w — kQW^p, (5-9) 
lyV^p^ lyk'^p- klp'^w. (5.10) 

They are the Euler equations determining the config- 
urations associated with the extrema of the Cole-Hopf 
Hamihonian i.e., 5H^^/5w = and 6H'^^/5p = 

0. By inspection we note that the Euler equations are 
compatible for p — (and w ^ 0) and for p oc w. For 
p = Eq. (|5.1U|) is satisfied identically, for p — vw Eqs. 
(|5.9(l and H5.10|l are identical; the prefactor v is dictated 
by dimensional arguments. 

On the noiseless manifold p = with H'^^ = we 
obtain the linear (Helmholtz-type) equation 

V^w = k^w. (5.11) 

An elementary solution of Eq. (|5.11|) is w oc exp(A:e • r), 
where e is a unit vector, = 1, pointing in an arbi- 
trary direction. This mode corresponds to the height 
field h — (k/koje ■ r, i.e., an inclined plane, and the con- 
stant slope field u — {k/kQ)e of magnitude k/ko point- 
ing in direction e. A general solution of Eq. (|5.11|) is 
constructed according to w{r) = J d'^e A{e) exp{ke ■ r), 
where since w > we must choose A{e) > 0. For the 
height field and slope field we obtain correspondingly 
h{r) — (l//co) log w(r) and u(r) — {k/ko)Vw{r)/w{r), 
respectively. By choosing the weight function A{e) ap- 
propriately we can prescribe the directional dependence 
of the exponential growth of w and the corresponding 
form of the height and slope fields. 

In the later analysis of the network solution it turns 
out that rotationally invariant solution of Eq. I|5.11|) . i.e., 



12 



an s-wave state, will be important in implementing the 
long distance boundary conditions. In polar coordinates, 
ignoring angular dependence, Eq. H5.11|l takes the form 



aw 



d — 1 dw 
r dr 



(5.12) 



At long distances, ignoring the first order term, we have 
w oc exp(±fcr). Incorporating the first order term by set- 
ting w ocr" exp(±fcr) and choosing the growing solution 
we obtain 

w° (r) - r^/^~d/2 exp(fcr), r ^ oo. (5.13) 

At small distances and in order to obtain a finite w at 
r = we must choose dw/dr — est. x r, implying 

wl{r) ~ w°_^{0) (X r'^ , r 0. (5.14) 

The exact solution of Eq. H5.12|) . finite at the origin, is 
given by w° (r) cx r-'^~'*/^/j;/2-i(fc''}jWhere Iu{z) is the 



Bessel function of the second kind 

Correspondingly, the asymptotic height and slope 
fields have the form 

hl{r)^^r, ulir)^~, r ^ oo, (5.15) 
fco ^ kor 

and at small distances 

/i°_(r)cxr2, u1(r)cxr, r ^ 0. (5.16) 

At large distances the height field forms a d-dimensional 
cone which at short distances becomes a d-dimensional 
paraboloid. The slope field has the form of an outward- 
pointing vector field of constant magnitude k/ko; for 
small r the slope field vanishes like r. In ID I_i/2{kx) = 
(2/7rA:a::)^/'^ cosh/cx and we obtain the Cole-Hopf, height, 
and slope fields in Eqs. (|5.6|l . (|5.5|l . and H5.4|) . i.e., the 
fields pertaining to the right hand domain wall solutions 
for the noiseless Burgers equation. 

On the noisy manifold p = vw we obtain the stationary 
nonlinear Schrodinger equation (NLSE) 



k^w — kiiif' , 



(5.17) 



which can be recognized as the stationary Gross- 
Pitaevski equation for a real Bose condensate with energy 
k^ and coupling strength fcg j^^- In radial coordinates 
we have 



aw 
dr^ 



d — 1 dw 
r dr 



^k\ 



(5.18) 



which supports a nodeless bound state w'^{r) falling off 
at large r and finite at the origin '98l|. For small w we 
recover the linear equation (|5.12() with decaying solution 



w°.{r) 



„l/2-d/2 



exp(— fcr). 



oo. 



(5.19) 



For small r, requiring finiteness of the first order term, 
we infer 



w°_{r) -w°_{0) (X -r^, r ^ 0. 



Correspondingly, the asymptotic height and slope fields 
have the form 

h°_(r) ^ --^r , u°(r)---^-, r ^ oo, (5.21) 
fco fco r 



and at small distances 



/i° (r) (X -r^ , u° (r) cx -r. 



0. 



(5.22) 



At large distances the height field forms an inverted 
d-dimensional cone, at short distances an inverted d- 
dimensional paraboloid. The slope field forms an inward- 
pointing vector field of constant magnitude k/ko; for 
small r the slope field vanishes like — r. In ID the radial 
equations H5.18|l takes the form d^w/dr^ — k'^w — k^w^ 
and admits the solution w^_{x) = (fc/fco)-\/2 cosh" kx in 
accordance with Eq. H5.5|l . yielding the height and slope 
fields in Eqs. (|5.5|) and (|5.4|l . i.e., the fields for the left- 
hand domain wall constituting the noise-induced growth 
modes for the ID noisy Burgers equation. 

By a simple scaling argument, w — > /xw, fco fco//^: 
the couphng strength fco can be scaled to 1. Conse- 
quently, the length scale is set by k~^ . Generally, 



uP_{r) = ad-^f{kr), 
fco 



(5.23) 



(5.20) 



where f{kr) ^ (fcr)*^^"'''/^ exp(— fcr) for kr '3> 1. Normal- 
izing f{kr) — 1 for r = 0, i.e., /(O) = 1, the dimensionless 
coefficient is a function of the spatial dimension d. In 
d = 1 we have from above ai ~ \/2, in higher dimension 
ad is determined numerically. We find 02 — 2.21 and 
03 = 4.34. In d > 4 the bound state is absent. This in- 
teresting feature will be discussed later in the context of 
the upper critical dimension. In Fig. |^ we have depicted 
the radially symmetric bound states of the NLSE for di- 
mensions d = 1, d = 2, c? = 3, and d = 3.5. The static 
localized spherical modes w° (r) and w^ (r) constitute the 
building block in establishing the growth morphology of 
the KPZ equation. In ID they become for the slope field 
the right hand and left hand domain wall solutions. In 
Fig. ^1 we have depicted the two kinds of radial modes 
for the w-field, /i-field, and u-field. 

Following Finkelstein ^| we finally give a simple argu- 
ment based on dynamical system theory for the existence 
a radially symmetric bound states of the NLSE 15.18|l . 
Treating r as a time variable Eq. 15.18|l describes the 
motion of a particle of unit mass with w as position in 
the double weU potential V{w) = ~{k^ /2)w^ + {k^/'i)w'^, 
subject to the time-dependent damping {d— 1) (dw/dr) /r. 
The phase space is spanned by w and dw/dr and char- 
acterized by a saddle point at {w, dw/dr) — (0,0) and 
an elliptic point at [w, dw/dr) = (fc/fco,0). The invari- 
ant homoclinic orbit (the separatrix) intersect the w axis 
at w = and w = ^/2k/ko. In Fig. 1111 we have in a 
phase space plot depicted the constant energy surfaces 
E = {l/2)(dw/dr)'^ + V{w) in the absence of damping. 
In d = 1 the damping is absent and the bound state solu- 
tion w-{x) oc cosh~^ kx corresponds to the motion along 
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FIG. 9: The radially symmetric bound states of the NLSE 
are depicted for k — ko — 1 in d = l,d = 2,d = 3,d = 3.5. In 
the limit d 4 the bound state solution disappears. 




> W_ 



FIG. 11: Phase space plot of dw-/dr vs. w_. We depict the 
constant energy surfaces E = ( 1/2) (dw/dr)^ + V{w) = est. 
The dashed line corresponds to the radial bound state solution 
of the NLSE. 




w_ 



■kr 




the invariant curve, the separatrix, from w — ^/2k/ko 
to the saddle point. For higher d, in the presence of 
damping, it follows from the equation of motion that 
dE/dr = — (d— l){dw/dr)'^ /r < 0, i.e., the energy de- 
creases monotonically in time. The orbit from an initial 
W-{0) with dw-/dr\r=o = will intersect the energy sur- 
faces and in most cases spiral into the elliptic fixed point 
at (fc/fco,0). Only a particular orbit, corresponding to 
the bound state, will reach the saddle point at (0, 0) for 
r — > cx). It follows from the numerical analysis that for 
increasing dimension, i.e., increasing damping, the initial 
value W- (0) defining the bound state solution migrates 
to larger values and out to infinity for d > 4. We have 
unsuccessfully attempted to determine the critical dimen- 
sion by a dynamical system theory argument; however, 
in Sec. IVIIII we determine the upper critical dimension 
algebraically by means of an application of Derrick's the- 
orem. 



\ 



2. Dynamical growth modes 




FIG. 10: We show the two static radial modes for the Cole- 
Hopf, KPZ, and Burgers fields. The KPZ fields correspond 
to cones at large distance and paraboloids at short distances; 
the Burgers fields are outward-pointing and inward-pointing 
vector fields vanishing linearly at the origin. 



The fundamental building blocks in establishing the 
growth morphology of the KPZ equation are the 
static localized modes Wj.(r), yielding the height fields 
h\{r) = {l/k[))\ogw^{r) and slope fields Uj_ — 
{l/kQ)^w±{r)/wj_(r). The growing mode is asso- 
ciated with the noiseless manifold p ~ and carries ac- 
cording to Eqs. (|4.9|) . (|4.13l) . and l|4.17|l no dynamics, 
i.e., = H'^^ ^ 0. The decaying bound state 

mode w'l, on the other hand, is endowed with dynamical 
attributes. The mode lives on the noisy manifold p — vw 
and is associated with the noise field p^_ = vvP, . Accord- 
ing to Eqs. (|4.13|) . H4.16|l . and 14.1711 it carries action. 
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energy, and momentum, 



CH 



H 



n^H = 0. 



(5.24) 

(5.25) 
(5.26) 

In order to generate dynamical modes we boost the static 
modes by means of the Galilei transformations H2.6|l , 
(I22|), l(TTT|l . and For the propagating locaHzed 

Cole-Hopf, KPZ, and Burgers modes we then obtain 

u;±(r, t) = w°^{\y + Au°i|) exp(fcouO • r), (5.27) 

(5.28) 

(5.29) 

The height field is a tilted d-dimensional cone mov- 
ing with velocity — Au*^; the slope field a d-dimensional 
hedgehog structure with an imposed constant drift 
u*^ moving with velocity — Au". In Fig. ^1 we 
have shown the height and slope fields in 2D. Note 
that in ID we have w^{x) oc cosh^ fcx, yielding 
hPj^(x) — (1/fco) log I cosh^ fexj -I- est. and Uj_(a;) = 
±(fc/fco) tanhfca;, and we obtain the propagating modes 



/i±(r,t) = 4(|r + Au"t|) + u°T, 
u±(r,i) =<(|r + Au"i|) + uO. 



discussed in 



u+{x,t) = ±(fc/A:o) tanh fc(a; 
Refs. ImIH 

The dynamics of the propagating modes is easily in- 
ferred. The noiseless mode (W+, u+) on the p — Q 
manifold has vanishing action, energy, and momentum. 
The noisy mode (w_,/i_,u_) on the p- = vw- mani- 
fold carries finite action, energy, and momentum. Since 
the action is invariant under the the Galilean boost the 
action is given by Eq. H5.24|l . For the energy and mo- 
mentum we obtain by insertion in Eqs. 1)4. 9|l and 14.17(1 . 
noting that the noise field is transformed according to 
Eq. (Ol^ . 



H 



CH 



1 



[kovf / d'^x w° (r)^ 



Un 



.0(r)^ 



n 



CH 



kQVUQ / d'^x w^_{rf 



(5.30) 
(5.31) 



Since the velocity of the mode is v = — Au° the ex- 
pressions (|5.30(l and H5.31|l admit a particle interpreta- 
tion. Defining the mass 

^ ^dW(r)2, 



we obtain 



H 



CH 



i?o^" -t- (l/2)m_^;^ 
where 

^CH ^ ^ii/2){kovf j d^x w°_{rY 
is the rest energy. 



(5.32) 

(5.33) 
(5.34) 

(5.35) 




FIG. 12: We depict the propagating height field and slope 
field modes in 2D. The bold arrows indicate the propagation 
velocity (arbitrary units). 



VI. 



PATTERN FORMATION 



The existence of the localized propagating growth 
modes in the Cole-Hopf field and the corresponding con- 
figurations in terms of the height and slope fields makes 
it natural to describe a growing interface in terms of a 
gas of growth modes. The implementation of this scheme 
constitutes a generalization of the weak noise approach 
in the ID Burgers case to higher dimensions. 



A. General properties 

From the general discussion in Sec. IIIII and Sec. IIVI it 
follows that a solution of the Cole-Hopf field equations 
()4.11(l and ((4.12(1 from an initial configuration ix;i(r) at 
time t = to a final configuration w{r) at time T, with 
the noise field p(r, t) as slaved variable, corresponds to a 
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specific kinetic pathway, yielding according to Eqs. (|4.13() 
and (|4.14|) the corresponding Arrhenius factor. The first 
central issue is thus to determine a global solution of the 
field equations. In th e spirit of instanton calculations in 
field theory (QOl llOCt HoiJ and following the scheme im- 
plemented in the case of the noisy ID Ginzburg-Landau 
equation and the noisy ID Burgers equation jOJ] we 
attempt to build a global solution on the basis of the 
propagating localized growth modes. In order to min- 
imize overlap contributions we, moreover, consider the 
case of a dilute gas of growth modes. In order to char- 
acterize a kinetic pathway or growing interface we must, 
moreover, choose appropriate boundary conditions. Here 
it is natural to assume a flat interface at infinity, that is a 
vanishing slope field. Note that this boundary condition 
does not preclude an offset in the height field and thus 
allows for a propagation of facets or textures accounting 
for the noncquilibrium growth. 



B. Dynamical pair mode 

In order to illustrate how to construct solutions we 
here consider a pair mode built from two growth modes. 
At time t = we combine a static noiseless and a static 
noisy Cole-Hopf mode centered at ri and r2, respectively, 
and obtain the three equivalent pair fields: 



Upair (r , 0) = (r - ri ) w° (r - r2 ) , 



(6.1) 



t-pair 



(r,0) = /i°(r-ri) + /i'L(r-r2), (6.2) 
Upai,.(r,0) = u0(r-ri) + u0_(r-r2). (6.3) 

For r — > oo we have u^, ~ {ki/ko)r/r and u*^ ~ 
— {k2/ko)r/r, where ki and k2 are the inverse wavenum- 
bers associated with the modes. Since Upair(r,0) ~ 
{{ki — k2)/ko)r/r we ensure a vanishing slope field at 
infinity by balancing ki and k2, i.e., ki — k2 ~ k. 

In order to assign velocities to the modes we note 
that in the vicinity of the mode position ri the slope 
field is shifted by u'?_(ri — Y2) and we must, accord- 
ing to the Galilei transformation (|2.11() . assign the ve- 
locity vi = — Au'^(ri — Y2). Likewise, the mode u'^ 
is assigned velocity V2 = — Au1(r2 — ri). For large 
separation |ri — r2| S> k^^ the asymptotic expressions 
(I5T5I) and i5^ yield Vi = \{k/ko){ri - r2)/|ri - Y2\ 
and V2 = —X{k/ko){r2 — ri)/|r2 — ri|- We note that 
Vi = V2 = V, i.e., the individual modes propagate col- 
lectively in a direction along the axis ri — r2 of the pair 
mode. The propagating pair mode is thus given by 

Wpah {r,t) = ^^(r - vt - ri)w" (r - vt - r2), (6.4) 
/ipair(r,0 = /i+(i--vt-ri) + /i"_(r-vt-r2), (6.5) 
Upair(r, t) = < (r - vt - n) + u'L(r - vt - ra), (6.6) 
v = -A<(r2-ri). (6.7) 

For the purpose of our discussion and assuming a dilute 
gas of growth modes it suffices to use the asymptotic 



forms of the Cole-Hopf, KPZ, and Burgers fields. Intro- 
ducing a core radius e of order 1/fc in order to regularize 
the solution at the origin and introducing the notation 



(6.8) 



we have for the propagating pair mode fields in more 
detail 

Wpair(r, i) exp(fc|r - ri{t)\^ - k\r ~ Y2{t%), (6.9) 



pa,r(r,t)~ ^(|r-ri(<)|,-|r-r2(t)|e), (6.10) 
fcn 



h 

Upair(r, t) 



"So 

k_ ( r-njt) 
ko ^|r-ri(t)|. 



r - r2{t) 
\r-r2it)U 



-X 



k r2 - ri 



ko |r2 - ril 

ri,2(0 =vt + Vx.2- 



(6.11) 

(6.12) 
(6.13) 



Denoting the separation d = r2 — ri and the unit vector 
f = r/r, and expanding h and u for large r we obtain 



k 

/lpair(l") = — f • d, 



u 



pair 



(r) = li[d-r(f -d)]. 
fco r 



(6.14) 
(6.15) 



We note that the height field is constant at infinity, its 
value depending on the direction, whereas the slope field 
vanishes. Introducing polar coordinates where ^ is the 
polar angle between the direction r and the axis d we 
have 



/ipair(r) = — dcos 
fc 1 

Upair(r) 



fco r 



d — fd cos ( 



(6.16) 
(6.17) 



In ID we obtain w\(x) cx cosh^^ fcx, /i±(a;) oc 
±(l/fco) log I coshfca;|, Wj.(x) = ±(fc/A:o) tanhfca;, and v = 
—XikjkQ) and we have for the pair mode the expression 

u{x, t) = 
k 

— [tanhA:(a; — xi{t)) — tanh/c(x — X2{t))], (6.18) 

fco 



-A- 



fco ' 

Xl.2{t) = Vt + Xi 



(6.19) 
(6.20) 



in accordance with Eq. H6.6|l and the discussion in Refs. 
[4^ : the ID pair mode is depicted in Fig. |S1 In Fig. 1131 
we have depicted the height and slope fields in 2D. 

The propagation of the pair mode is in the direction 
of the axis connecting the two centers ri and ri . During 
the passage of a pair with velocity v — Afc/fco it follows 
from Eqs. (|6.16() and H6.17|l that the local height and 
slope fields change by Sh = 2fcc?/fco and 6u — 2fc/fco, 
respectively. The passage time of the pair is St = d/v = 
dko/Xk and we obtain for the growth velocity Sh/St = 
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FIG. 13: We show the height and slope field of a dipole or 
pair mode in 2D. The bold arrows indicate the propagation 
velocity (arbitrary units). 



2A(fc/A:o)^ = in accordance with the averaged 

KPZ equation in a stationary state, {5h/5t) = (A/2)u^. 

The dynamics of the pair mode is inferred from Eqs. 
(|5^ . lfOn|l . and ifOT]! . Inserting Eq. lf05|l a scahng 
argument yields 



1 



d^S, (6.21) 



for the action of a pair propagating in time T. Note 
that only the decaying bound state component carries 
action. In the dynamical phase space language this sce- 
nario corresponds to a pair-orbit from an initial con- 
figuration ?i;i(r,0) — uP^{y — Yi)uP_{y — Y2) with initial 
noise field pi(r, 0) = vw^_{y — Y2) to a final configura- 
tion Wf(r,T) = uP^{y — vT — ri)u>" (r — vT — Y2) and 
p{{r,T) — i/it;° (r — vT — r2). In Fig. El we have in a 
{WjP) phase space plot sketched the specific orbit from 
(wijPi) to (wfjPf). Choosing the centers ri and r2 and 
the amplitude k we thus have a whole class of orbits cor- 
responding to kinetic transitions from (ri,r2,A:) at time 
i = to (ri -I- vT, r2 -t- vT, k) at time T; note that the 



FIG. 14; Phase space plot of a dipole or pair orbit from an 
initial configuration Wi to a final configuration W{ in transition 
time T. The noise field p is a slaved variable going from pi to 

PC- 



magnitude of the velocity is t; = Xk/ko, whereas its di- 
rection is given by d = ri — r2. 

For the rest energy and mass of the pair mode we have, 
inserting Eq. (jF?^ in Eqs. (|F3n|) . fOT| . fO^ . and 
(ESS), 



E, 



'0,pair 



C /(0',(6.22) 



1 / k 



and we obtain for the energy and momentum 



1 



pair I ^^pair"^ 



n 



pair 



''pair * 1 



ko 



(6.23) 



(6.24) 
(6.25) 
(6.26) 



The pair mode satisfying the boundary conditions of 
an asymptotically flat interface suggests an independent 
particle picture of the growth morphology in terms of a 
dilute gas of pair modes. The pair modes have masses 
scaling with the amplitude according to Wpair « fc^^'*. In 
2D the mass is independent of k, in ID the mass grows 
linearly with fc, for d > 2 the mass vanishes for large k. 
The rest energy i?o,pair oc fc'*^'*, i.e., i?o,pair is indepen- 
dent of A: for d = 4. For d < 4, i?o,pair grows with k, for 
d > 4, £^o,pair vanishes for large k. Finally, the velocity 
scales linearly with k. 



C. Dynamical network 

The generalization to a network of propagating modes 
is straightforward. At time i = we assign a set of 
growing and decaying static modes at the positions r?, 
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i = 1,- ■ ■, and obtain the Cole-Hopf field at time t = 0, 
w{r,Q) = l[wUr-r",). (6.27) 

i 

It is here convenient to use a charge language for the 
amplitudes ki pertaining to the i-th mode. For ki > 0, a 
positive charge, w,^(r) is a growing solution of the linear 
equation H5.11|l : for ki < 0, a negative charge, u'-'(r) is a 
decaying bound state solution of the NLSE H5.17|l . For 
large r we have w°{r — r-*) ^ exp(fcir), i.e., w{r,0) ~ 
exp(r^^fci) and in order to ensure an asymptotically 
flat interface h{r) — (l/ko) logw(r) we impose the charge 
neutrality condition 



^ fc, = 0. 



(6.28) 



The corresponding initial height and slope field are then 
given by 



h{r,0)^^^Y.^ogw°{r-v'',), (6.29) 

(6.30) 



ko 
1 

ko 



u(r o) = lyZ!fcii) 



In order to assign velocities to the modes we proceed 
as in the case of the pair mode. At the position r° of 
the £-th mode the slope field is shifted by u(r",0) = 
(1/fco) J2i Vw°(r^ - r?)/w;°(ry - r^), corresponding to a 
Gahlei boost to velocity = Au(r°, 0). Since, unhke the 
case of a pair, the modes move relative to one another 
the network will converge towards a self consistent state. 
We thus obtain the self consistent dynamical network 



w{r,t) ^l[w°{r-r,{t)), 



nit) = J w,{t')dt' + vl 
E fc, = 0. 



(6.31) 

(6.32) 
(6.33) 
(6.34) 



The interpretation of the growth morphology or pattern 
formation represented by Eq. lt).31l().34|) is straightfor- 
ward. In the weak noise WKB representation the grow- 
ing interface is described by a gas of growth modes with 
negative and positive charges. The asymptotic flatness 
condition is ensured by imposing charge neutrality as ex- 
pressed by Eq. H6.34|l . The dynamics of the network is 
constrained by the assignment of velocities to the modes, 
where according to Eq. Ht).33|l the velocity of a particular 
mode depends on the position and charges of the other 
modes. The connectivity and continuity of the network 
thus defines the temporal evolution. 

For the present purposes it is sufficient to consider a 
dilute network and use the asymptotic form u'^(r) ^ 



exp(fcir). We then obtain the height field, slope field, 
and assigned velocities 



n{t)-vt{t) 
\n{t)-v,{t)l 



(6.35) 
(6.36) 
(6.37) 



where we have used the short distance or UV regulariza- 
tion given by Eq. Ht).8|l . 

Since we initialize the network at rest, it will pass 
thorough a transient period where the velocities adjust 
to constant values as the modes recedes from one an- 
other. From Eq. H().37|) we obtain by differentiation 
dvi/dt = -2i^X]£5ziJ(vi - V£)/|rj - ti\^ - (r* -r^) • (v^ - 
we){vi ~ re)/\r,, - Yi\l] which, assuming - ^f^ to be 
bounded, vanishes for jr^ — r^j oo. We infer that at 
intermediate times longer than the transient time the ve- 
locities attain constant values. Using Yi{t) = v^t -I- we 
thus obtain from Eq. (|6.37|l a self consistent equation for 
the velocities in the stationary state 



-2v 



At fixed time for r large expanding Eq. 
tain 



h{Y,t) 



1 r 

fco r 



(6.38) 



we ob- 



(6.39) 



where we have used the neutrality condition (|().34|l . i.e., 
an asymptotically flat interface. 'Y^^ kiYi{t) defines a 
center of mass position R(t) = J^i^i^iW ^^r a di- 
lute cluster of modes. Introducing the polar angle (f> 
between the direction r/r and R(t) we have, in anal- 
ogy to Eqs. H6.16|l and (|6.17|) in the case of the pair 
mode, h{r,t) — — (l/fco)|R(i)| cos0, and the height off- 
set depends on the direction. As the cluster of modes 
propagate across the system the height field changes by 
2|R(t)|/fco. Likewise, for the slope field u expanding Eq. 
(|6.36() for large r we have u(r,t) ~ (r/?") X^i ''^j — 0' 
i.e., asymptotically a fiat interface. We also note that 
K{t) = j:^hir° + v,<) = E.fc.(r° + v.i) = E^^rl 
where we have used J2i^i'^i — 0, following from Eq. 
(|6.38l) . i.e., R(i) — J2i ^i^i independent of time in the 
stationary state. Likewise, at fixed position for t large we 
have by expanding Eq. Ht).35|l . h{r,t) — [t/koj^^kiVi, 
and we infer the constant growth velocity 



dh 



ko 



(6.40) 



The relationship between the imposed drift F in the 
KPZ equation (|l.l|l and the charge assignment to the 
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growth modes, {ki}, is given by 

^-^(Efc^) ■ (6-41) 

In order to demonstrate this identity we consider the 
Cole-Hopf field equation (|4.28|l in the asymptotic re- 
gions where the noise field vanishes, dw/dt = i^V^w — 
vk'^w, fc^ — \F/2v^. From the growth mode ansatz 
w — Hiexpfcilr — Yi{t)\ we readily obtain {dw/dt)/w — 
-Y^tki^^ ■ (r - i"j(t))/|r - r,(t)| and V'^w/w = 
Ei,- fc, (r - r, ) ) • (r - r, (i ) ) / 1 r - r, ) 1 1 r - r, (t ) I . Insert- 
ing the mode velocity from Eq. (|6.37() and symmetrizing 
we have for large r the identity (|6.41|) . The drift F is thus 
given by the total charge of the growth modes. A neu- 
tral system with vanishing slope at infinity corresponds 
to F = 0. 

In addition to the dynamical velocity constraint im- 
posed by the continuity and connectivity of the network, 
the canonical structure of the WKB weak noise scheme 
impart dynamical attributes to the network. As far as the 
dynamics is concerned only modes with negative charge 
on the noise manifolds pi = I'Wi, corresponding to the 
bound state solution of the NLSE, contribute. Accord- 
ing to Eq. (|5.24() we obtain for the total action of the 
network 

Snot = J2 (•^•42) 

ki<0 

where the action of the i-th mode is 




(6.43) 



For a dilute network using the asymptotic expression 
(|6.21l) we obtain accordingly 

5. = ^ikoi^fT (^^y ajk-^ J^°° d'C fiO'- (6.44) 

At time i = we assign an initial state Wi{r,0) by 
choosing a set of positions and associated charges ki . 
The expressions (|6.37ll and Ht).32l) subsequently determine 
the appropriate propagation velocities and time depen- 
dent positions. The network develops dynamically defin- 
ing a particular kinetic pathway. In order to obtain a 
specific path from an initial state wi to a final state W{ 
traversed in time T we must choose a specific set of po- 
sitions and charges. The action associated with the path 
is then given by Eqs. (|6.42I6.44|I . yielding the Arrhenius 
factor P oc exp[— S'/A] for the corresponding transition 
probability. 

Associated with the initial configuration wi at time 
t = is the initial noise field 

Pi(r,0)cx n ^"(i--r°). (6.45) 

ki<0 




W 



FIG. 15: Phase space plot of a set of orbits with the same 
transition time T terminating in different final configurations 

Assigning velocities the noise field propagates according 
to 

p(r,i)cx n ^i'°(i--v.t-rO), (6.46) 

ki<Q 

to the final noise configuration at time T. For a fixed 
transition time T and a given initial lu-configuration dif- 
ferent assignments of ki and yield different noise field 
and, correspondingly, different final {w,p) configurations. 
The scenario is depicted schematically in Fig. [T51 

The energy and momentum of a network is inferred 
from Eqs. H6.22I6.26|) 

End = Enet,0 + ^ X! "^^^^ (6.47) 

i 

n„ct = ^ rriiVi, (6.50) 

ki<0 

where in the stationary state is given by Eq. H6.38(l . 

D. Field theory 

In a qualitative sense the Cole-Hopf field w and the 
associated noise field p governed by Eqs. H4.11|) and 
(|4.12ll representing the KPZ equation in the weak noise 
limit plays the role of bare fields, whereas the propa- 
gating noiseless and noisy modes Wi together with the 
associated noise fields connected in a dynamical network 
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constitute the renormalized fields. The strong coupling 
features of the problem are thus represented by the net- 
work. It is instructive to represent this insight in terms 
of a field theory for the network. Below we sketch aspects 
of such a field theory. 

We consider a dilute distribution of modes or 
monopoles with charges ki at positions and introduce 
the density n and charge density p according to 

n(r,t)=^,5'^(r-r,(t)), (6.51) 

i 

p{v,t)^Y.^y{y^v,{t)); (6.52) 

i 

we note that the neutrality condition fcj = is equiv- 
alent to / d'^xp{v, t) — 0. For the asymptotic height and 
slope fields \iS.'6b\ and (|6.36l) we obtain 

h{v, t)^^ [ p{T',t)\v - r'|„ (6.53) 

u(r, t) = ^ / p(r', t)j^—^. (6.54) 



Introducing the velocity field 

v(r,t)=^v,(0<5'^(r-r,(t)), 



The equation of motion for the charges is given by Eq. 
(|6.56l) which determine the velocity field v in terms of 
the slope field u. It is instructive to compare Eq. (|6.55|l 
to the Lorentz equation for the motion of charges e with 
density n in an electric field E, mdv/dt — enE. Including 
a damping term —7V and considering the overdamped 
case we obtain v = — (e/7)nE which has the same form 
as Eq. (|6.56|) with the slope field playing the role of an 
electric field. 

Finally, introducing the continuity equations express- 
ing number and charge conservation 



dn 

'at 

dp 
dt 



- V • V = 0, 
V-j-0, 



where the current density is 



j(r,i) =^fc,v,,5'^(r-r,(t)), 



(6.62) 
(6.63) 

(6.64) 



(6.55) 



we have the field equation H6.61|l for u and the Lorentz 
equation (|6.56() for the particle dynamics. In Fig. 1161 we 
depict a neutral three-mode configuration in the height 
field and the associated divergence of the slope field. V-u 
clearly shows the two negative charges and the single 
positive charge constituting the morphology. 



we express the velocity condition (|6.37|) in the form 

v(r, t) ^ -2vn{Y, i)u(r, t). (6.56) 

Using V • (r/r) — [d — \)/r and introducing the scalar 
field or potential (j) 



= V • u = V^/i, 
we have from Eq. (|6.55() 



fco 



r - r' 



(6.57) 



(6.58) 



i.e., (j) is analogous to the electric potential arising from 
a charge distribution p{v,t) [loij . It also follow from Eq 
(|6.58|) that 4> satisfies the fractional Poisson equation 



V 



d-i 



d-1 
— / 



(6.59) 



where V''^^ is the Fourier transform of k'^^^; note that 
in d = 3 Eq. H6.59II becomes the usual Poisson equation 
in electrostatics 102]. 

The charge distribution p yields the potential <j) either 
as a solution of Eq. H6.59|l or in terms of the integrated 
form in Eq. H6.58|l . Inserting Eq. H6.57|) we also have 



V^-^V -ucx 



d-l 



ko 



-P, 



d~l 
ko ' 



(6.60) 
(6.61) 



E. Linear diffusive modes 

Generalizing the discussion in the case of the ID noisy 
Burgers equation to higher dimension, it is clear that in 
addition to the network of growth modes there are also 
superimposed linear diffusive modes. Here we summarize 
aspects of the linear mode spectrum. 



1. The linear case 

In the linear Edwards- Wilkinson case the weak noise 
field equations H4.20|l and H4.21|l take the form 

//II 

|^ = i.V2u-V(V.p), (6.65) 

^ = -i/V^u. (6.66) 

Since u(= Vft-) is longitudinal and since only the longitu- 
dinal component of p couple to u, Eqs. (|6.65(l and (|6.66|l 
for the wavenumber components Uk and pk correspond 
to the overdamped oscillator case discussed in appendix 
IXI For an orbit from uj^. at time t = to uj^. at i = T, 
with Pk as slaved variable, we have 



Uk(t) 

Pk(t) 



Uj^ sinh LJk^ + sinh cJk(2^ — t) 



ve 



sinh LO\^T 



sinhcjk?^ 



, (6.67) 
(6.68) 
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FIG. 16: We depict a neutral three-mode configuration in the 
height field and the associated divergence of the slope field 
(arbitrary units). 



with the diffusive mode frequency 



(6.69) 



The spectrum is exhausted by linear diffusive modes with 
eapless dispersion given by Eq. (|6.69|) . Following Ref. 
|42| the action and transition probabilities are given by 



d'^k |u(-uj^exp(-u;kr)|^ 



(27r)'' 
P(u(,uL,T)aexp 



1 - exp(-2wkT) 
S' 



(6.70) 
(6.71) 



In the limit T ^ oo we obtain the stationary distribution 

dfik . 



Po(uk) oc exp 



'a 



Uk 



(6.72) 



yielding the Boltzmann distribution H2.1|l : for more com- 
ments on the linear case see the discussion in Sec. m 



2. The nonlinear growth case 

Iir the nonlinear KPZ case the growth morphology or 
pattern formation is given by a dynamical network of 
propagating localized growth modes. Superimposed on 
the dynamical network is a spectrum of extended linear 
diffusive modes. In order to implement the boundary 
condition of an asymptotically flat interface it is most 
convenient to conduct the discussion in terms of the slope 
field u. 

Considering the field equations (|4.20(l and H4.21|l we 
set u — Unot + <5u and p = p„ct + Sp. Here Unet is given 
by iniSnil, i-e., 



u„ot(r,t) = T- y]fc: 

Kg 



(6.73) 



For the associated noise field Pnot we obtain from 
Eqs. H4.33|l . H5.19|l . and (|4.34() for a single mode p = 
kopw = k^vuPi cx k[)vr^^'^ eyiY){—2kr), i.e., V • p = 
— fcoi^r^"'^ exp(— 2A:r). Setting p = {Y/r)g{r) we have 
dg/dr + ((d— l)/r)g — —k^vr^^'^ exp(— 2fcr) with asymp- 
totic solution g{r) ^ {koi^ /2k)r^^'^ exp{~2kr). General- 
izing to a propagating network we have 



Pnet(r,i) 



_iykos-^ r-r,{t) 2|fc.||r-r 



fei<0 



l^.l|r-r,(i)|^ 



(t)l 



(6.74) 



We note that the noise field is localized in the vicinity of 
the modes with negative charge, i.e., the modes carrying 
dynamics. The noise field falls off exponentially with 
a range given by Confining our discussion to the 

regions between the growth modes and noting that Pnot 
and Vunct ~ we obtain linear equations for i5u and 
Sp, 



dSu 
dSp 



= iyV^6u + A(u„ct • V)(5u + V(V • dp), (6.75) 



-J/V^Jp -I- A(Unct • V)6p. 



(6.76) 



Assuming that u„ot ~ const, and setting 5u, 5p ^ 
exp{—Et) exp(ikr) it follows that 

Su = (ae-f'^'^'-'^"-'-'')* + be('''='+*^""■='■'')*)e*''^(6.77) 
(5p = ce('^'='+*^"-'•'')*e'''^ (6.78) 

For A = we recover the linear case; in the nonlinear 
case, incorporating the time dependent term iAu,jot ■ kt 
in the plane wave phase ikr, we note that the diffusive 
modes undergoes a mode transmutation to damped and 
growing propagating modes with phase velocity Aunot • k. 



VII. 



STOCHASTIC INTERPRETATION 



The weak noise scheme gives access to the transition 
probability P{hi /if, T) from an initial height profile 



21 



hi to a final profile /if in time T, where P is given in 
terms of the action S{hi h{,T), P cx exp[— 5'/A]. We 
note that the scheme only yields what corresponds to the 
Arrhenius factor exp[-S'/A]; the prefactor r(r), P{T) = 
r(r) exp[— S'(T)/A], is to leading order in A determined 
by the normalization condition J J^^. dh{P{hi —> /if, T) = 
1, i.e^r-i(T) = /rir'^^f^eii huT), see also Refs. 
M 



A. Kinetic transitions 

The prescription is in principle straightforward. An 
initial configuration hi at time t = is modelled or ap- 
proximated by a dilute gas of growth modes forming the 
network characterized by their charges {ki\ and posi- 
tions {r°} plus a spectrum of diffusive modes {5u\,5p\}. 
The dynamical configuration evolves in time according to 
the field equations. The velocities of the growth modes 
are assigned according to Eqs. H6.33|l and H6.37|l : notice 
that the noise field associated with the negatively charged 
modes develops in time according to Eq. H6.74|l . At time 
t = T the profile h has evolved to the final profile /if, cor- 
responding to the network configuration {ki,Ti{T)} and 
the final diffusive mode configuration {^u^,<5p^}. This 
time evolution corresponds to a specific kinetic pathway 
for a growing interface. The action associated with the 
transition is composed of a network part Snet and a dif- 
fusive mode part S'difT- Snct is given by Eqs. Ht).42|l . 
(|6.43|) . and Ht).44|l and thus only depends on {ki} for the 
negatively charged bound states or monopoles. In the 
linear case S'difT is given by Eq. (|6.70|l : we note that 
SdiS depends on the initial and final diffusive mode am- 
plitudes 6u\. and Su^,. At long times Snot grows linearly 
with T, like in the case of random walk discussed in ap- 
pendix ^ whereas Sdiff approaches the stationary form 
L' J d'^k/{2n f\ul.\'^, see Eq. (|^7^ . as in the case of the 
overdamped oscillator discussed in appendix IXI 



B. Anomalous diffusion and scaling in the dipole 
sector 



Leaving aside the issue of the linear diffusive modes, 
the network representation is based on an assumption of a 
dilute gas of growth modes or charged monopoles. In the 
course of time the modes will in general collide and coa- 
lesce and the dilute gas approximation ceases to be valid. 
Here we consider the class of network configurations com- 
posed of a dilute gas of pair modes or dipoles. A single 
dipole satisfies the boundary condition of vanishing slope 
field. Consequently, the dipoles move independently. 

A single dipole or pair mode with charges k and ~k 
propagate according to Eq. H6.26|l with velocity 



k 

ko ' 



(7.1) 



and carries according to Eq. H6.21|l the action 
S = Tk^-d( 

ko 



(7.2) 



where A = d'^S, /(^) only depends on dimension. 

During time T the center of mass of the dipole propagates 
the distance L = vT and eliminating the charge k we 
obtain the action 



(7.3) 



where B = 

The form of Eq. 1)7.3(1 allows an interpretation of 
the ballistic motion of the dipoles within the dynamical 
scheme as a random motion within the stochastic descrip- 
tion. From the WKB ansatz we obtain for the transition 
probability P(L,T) over a distance L in time T for a 
single dipole 



P(L, T) cx exp 



B L 



4-d 



A ps-d 



(7.4) 



By a simple scaling argument the mean square displace- 
ment is given by 



T 



2H 



2{1-H) 



T{3-2H), 



(7.5) 



where T{z) is the Gamma fu nction |96 
troduced the Hurst exponent |l03j 



H = 



3-d 



and we have in- 



(7.6) 



From the scaling form in Eq. (|2.13|) and from Eq. H7.4|l 
it also follows that the dynamic exponent for the dipole 
sector is z = ^/H, i.e.. 



4-d 
3-d' 



(7.7) 



In ID we obtain H = 2/3 and 2: = 3/2 in accordance 
with well-established results for the KPZ equation or, 
equivalently, the noisy Burgers equation. Here the dy- 
namic exponent for the dipole sector agrees with the ex- 
act value. Note that the scaling law H2.14|l . z + ( — 2, 
following from Galilean invariance and automatically im- 
plemented in the present approach, implies the roughness 
exponent ( = 1/2. The Hurst exponent H — 2/3 > 1/2 
corresponds to a persistent random walk. Since the mean 
square displacement for i7 = 2/3 > 1/2 falls off faster 
than Brownian walk (H = 1/2) we have the case of su- 
perdiffusion. In 2D, which in a scaling context is the 
lower critical dimension, wc have H = 1/2, corresponding 
to ordinary Brownian diffusion. The dynamic exponent 
z = 2 and according to z + C = 2 the roughness exponent 
C = 0, corresponding to a smooth interface. In 3D the 
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VIII. UPPER CRITICAL DIMENSION 



1 



1/2 



H 
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- diffusion 



a) 



sub- 
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b) 



Brownian 
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log 



logT 



frozen 



FIG. 17: In a) we plot the Hurst exponent as function of d. 
In b) we depict in a log-log plot the mean square deviation 
(51/^) as a function of T. 



Hurst exponent H = corresponding to the logarithmic 
case 



{6L^) cx logT for d = 3. 



(7.8) 



Here the mean square displacement falls off slower than 
Brownian diffusion and the random motion of the dipole 
modes is characterized by being antipersistent and show- 
ing subdiffusion. For 3 < d < 4 the Hurst exponent 
H < and the mean square displacement decays cor- 
responding to pinning or glassy behavior. In 4D, which 
we here propose to be the upper critical dimension for 
the scaling properties of the KPZ equation, the Hurst 
exponent H = — oo, corresponding to extreme pinning or 
arrested growth. In Fig. 117b we depict the Hurst expo- 
nent as a function of dimension, in Fig. 117b we plot the 
dipole mean square displacement as function of time in 
a log-log plot. 



In addition to the scaling properties in the rough phase 
characterized by a strong coupling fixed point which in 
ID yields z = 3/2, a major open problem is the existence 
of an upper critical dimension. On the basis of the singu- 
lar behavior of perturbation theor y an d the beta function 
in a Callan-Symanzik RG scheme j4J|, the mapping to di- 
rected polymers 45] , and mode coupling arguments , 
it has been conjectured that d = 4 is an upper critical 
dimension for the KPZ equation. The behavior above 4D 
presumed to be complex and maybe glassy is, however, 
not well understood. 

Here we address the issue of an upper critical dimen- 
sion within the context of the weak noise approach and 
associate it with the existence of growth modes. This 
is not a scaling argument but based on the assumption 
that the growth mechanism and strong coupling features 
of the KPZ equation depend on the existence of prop- 
agating localized growth modes across the system. In 
Sec. the numerical analysis of the bound state solu- 
tion of the NLSE indicated that for d > 4 the solution 
is absent. This suggest that d = 4 within the present 
interpretation plays the role of an upper critical dimen- 
sion. Below we amplify this argument by means of an 
application of Derrick's theorem. 



A. Derricks theorem 

Derrick's theorem 0, see also Refs. [Toilliol . 
states that for a wide class of nonlinear wave equation 
there does not exist stable localized time- independent so- 
lutions with finite energy in dimensions greater than one. 
The theorem effective rules out soliton-like solutions to 
Lagrangian field theories in higher dimensions. Here we 
sketch the simple arguments in Derrick's theorem. 

Let us consider the generic Hamiltonian 

H=^l d''x[{Wcl>f+V{cb)+p% (8.1) 

for a scalar field in d dimensions. V((j)) is the potential 
and p the momentum. From the Hamilton equations of 
motion d(t>/dt = 5H/5p and dp/dt = —5H/S(j) we readily 
obtain d(j)/dt = p and dp/dt = V'^cj) - (l/2)y'(0), or 
eliminating p, the nonlinear wave equation 



d 



(8.2) 



For a stationary field 0(r) the energy is given by 

E=\jd''x[[V<l>f + Vm, (8.3) 
and the Euler equation 

V20=ir(0), (8.4) 
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follows from the variational principle 
5E 



(r) 



= 0. 



(8.5) 



To ensure stability we, moreover, require that the matrix 



5^E 



1 



5(t){v)5(t>{v') 



V2 + -T/"(0)](5(r-r'), (8.6) 



is positive definite, i.e., S^E > 0. Performing a con- 
strained minimization corresponding to a dilatation or 
scale transformation 



and introducing the notation 



(8.7) 



I = \j d'^xVicj,), (8.9) 
E^ = \j d^xliVcP^f + V{c^^)l (8.10) 



we obtain by substitution the relation 



Ef, = ^ ^-aj 



-d-, 



.11) 



Implementing the variational principle SE/dcj) = as a 
constrained variation dE^/dii\fj,^i = we infer the iden- 
tity 



(2 - d)K = I. 



.12) 



Moreover, the stabihty matrix 6'^E/S(j)'^ ~* 
is given by 



d^E, 



—f] ={2-d}{l-d)K + d{d+l)I. (8.13) 



Inserting Eq. H8.12|l we finally obtain 

'd^E. 



^ ' = 2(2 - d)K. 



.14) 



Since if > we must require d < 2 in order to obtain a 
stable stationary solution. 

This is the basic result of Derrick's theorem: The ex- 
istence of stable localized configurations (/)(r), i.e., with 
a finite norm, as solutions to nonlinear field equations, 
is only ensured below 2D. Derrick's theorem is a no-go 
theorem which effectively rules out stable soliton-like so- 
lutions in higher dimension. 



B. Nonexistence of growth modes above d = 4 

Here we address the issue of the existence of growth 
modes by an application of Derrick's theorem to the 



NLSE, see a lso Refs. [106|, [lO'Tl- In the context of the 
NLSE (|5.17|) the issue is not stability but existence of 
bound states. Here the energy functional can be ex- 
pressed in the form 



E = K +-k^N -'^L 
2 4 



.15) 



where K is the bending energy, N the norm, and / the 
interaction, 

K=\ j d'^xiVwf, (8.16) 

(8.17) 

/= / d'^xw^, (8.18) 



N = I d'^xw^ 



and the NLSE 



(8.19) 



follows from the variation 6E/6w — 0. 

In order to establish the bound we need two identities. 
The first identity is obtained by multiplying the NLSE 
by w and integrating over space yielding 



2K + k^N-kll^0. 



.20) 



The second identity is obtained by constrained minimiza- 
tion. Subject to the scale transformation w{r) — ?x;(/ir) 
we infer K /^^"^K, N ^'^N, and / /^^/, and ap- 
plying constrained minimization dE/dfi\fj_^i — we infer 
the second identity 



k^ k^ 
(d - 2)K + —dN ~ ^dl = 0. 
2 4 



.21) 



Eliminating the bending term K from the two identities 
we obtain 



fc2Ar^!J(4-d)/. 



.22) 



Since > 0, and / > it follows that d < 4 in order for 
a bound state to exist. 

This completes the proof of the nonexistence of bound 
states and thus growth modes for the KPZ equation 
in the Cole-Hopf formulation in dimensions larger than 
4. The proof corroborates the numerical analysis of the 
bound state solution. 



IX. SUMMARY AND CONCLUSION 

In the present paper we have extended the weak noise 
approach, previously applied in detail to the ID noisy 
Burgers equation, to the KPZ equation in higher dimen- 
sions. Three issues have been addressed: i) kinetic pat- 
tern formation, ii) anomalous diffusion and scaling, and 
iii) the upper critical dimension. 
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i) The weak noise WKB formulation allows a classi- 
cal interpretation of the pattern formation in the KPZ 
equation in the sense that the growth of the interface 
is interpreted as a dynamical deterministic network of 
propagating localized growth modes. The growth modes 
play the role of elementary excitations and are analogous 
to the vortex structures in the Kosterlitz-Thouless the- 
ory or hedgehog structures in the ferromagnet, see Ref. 
0. The imposed network structure expresses the strong 
coupling features. Superimposed on the network is a gas 
of subdominant extended diffusive modes corresponding 
to the EW universality class. The dynamical evolution of 
the network together with the diffusive modes defines the 
kinetic pathways from an initial configuration to a final 
configuration. Within the canonical weak noise scheme 
the network is endowed with dynamical attributes, it car- 
ries energy, momentum, and action. Here the action S 
plays the particular role of a weight function in deter- 
mining the transition probability P cx exp(— 5/A), A is 
the noise strength, for a specific kinetic pathway. 

ii) The weak noise method gives access to the scal- 
ing properties of the KPZ equation. The nonperturba- 
tive character of the WKB approximation implies that 
strong coupling features might be accessible. This is in 
fact the case in ID where the exact scaling exponents 
can be retrieved from the dispersion law for the growth 
modes and the structure of the stationary zero-energy 
submanifold. In higher D the stationary submanifold is 
not know and only limited scaling results are available in 
the dipole sector. In the dipole or pair mode sector, cor- 
responding to a dilute gas of dipole modes, the stochastic 
interpretation implies that the pair modes perform ran- 
dom walk with Hurst exponent H — (3 — d)/ (4 — c?) . In 
ID the interface grows stochastically with H = 2/3, cor- 
responding to persistent super diffusion; in 2D, the lower 
critical dimension, H = 1/2, corresponding to ordinary 
Brownian motion; in 3D we have H = 0, equivalent to 
logarithmic antipersistent subdiffusion; finally, in 4D, the 
conjectured upper critical dimension, H diverges, corre- 
sponding to an arrested or frozen interface. Formally, 
the dynamic exponent for the dipole sector is given by 
z = 1/H = (4 — c?)/(3 — (i). In ID we recover the well- 
known result z = 3/2; in 2D we have z = 2 which is 
the weak coupling result. In 3D the dynamic exponent 
diverges; in 4D we have z — 0. We believe that the be- 
havior of z above c? = 2 is an artifact associated with the 
dipole sector; it does not reflect the true scaling behavior 
of the KPZ problem. 

iii) The issue of an upper critical dimension for the 
KPZ equation has been much debated. Dynamical renor- 
malization group arguments indicates d = 4 as the upper 
critical dimension; numerical simulations, on the other 
hand, suggest an infinite upper critical dimension. Here 
we associate the upper critical dimension with the exis- 
tence of growth modes. This is not a scaling argument 
but rather associated with the idea that in the absence 
of a network of localized growth modes there is no clear 
mechanism for the growth of the interface. Above 4D 



the interface is, of course, still governed by the nonlinear 
KPZ equation but we lack insight into the actual growth 
mechanism if any. 

This paper constitutes our present understanding of 
the KPZ equation within the weak noise WKB approach. 
However, many open problems remain. 

i) Dynamical network: We have only analyzed the net- 
work in terms of the asymptotic form of the growth 
modes. A detailed analysis of the long time behavior 
requires further analysis of the field equations in order to 
understand collisions and coalescence of growth modes. 
Moreover, the analysis of the linear diffusive modes and 
their interaction with the growth modes is incomplete. 
Further analysis is required in order to determine the 
transition probabilities for specific kinetic pathways. 

ii) Scaling: In order to determine the scaling properties 
we need to determine the stationary zero-energy manifold 
yielding the stationary distribution. This is only possible 
in ID; in higher D the stationary submanifold is unknown 
and we only have scaling results for the dipole sector. 

iii) Dynamic renormalization group: The present ap- 
proach is based on the noise strength being the nonper- 
turbative small parameter in the problem. In that sense 
the method is not a scaling approach. It is unclear how 
to relate the weak noise approach to the dynamic renor- 
malization group scheme. 

iv) Directed polymers: An important element in our 
understanding of the KPZ equation is the mapping to 
directed polymers in a random medium. The relation- 
ship between the growth modes and the wandering of 
polymers is an open problem. 

In conclusion, the present nonperturbative weak noise 
scheme represents an alternative angle of approach to the 
KPZ equation and similar noise-driven problems. We em- 
phasize that the method is not a scaling approach based 
on expansions about critical dimensions, but rather iden- 
tifies the noise strength as the small parameter in the 
problem. Moreover, the method is not based on per- 
turbation theory but on an asymptotic nonperturbative 
WKB or eikonal approximation. 
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APPENDIX A: WEAK NOISE APPROACH TO 
SIMPLE RANDOM WALK AND THE 
OVERDAMPED OSCILLATOR 

1. Random Walk 

Simple ID random walk is described by the Langevin 
equation 

^^vit), {vit)m) = m), (Ai) 

corresponding to vanishing drift F ~ 0. From the weak 
noise scheme we obtain the Hamiltonian 



H 



1 



and the equations of motion 

dx dp 

with solution from to x in time T, 

x-x" ^Pf)T , p= Po- 



rn 



(A3) 



(A4) 



In the weak noise scheme random walk corresponds to 
free particle propagation. The action is 



and we infer the well-known distribution 

(x - x' f 



P{x, T) cx exp 



2AT 



(A5) 



(A6) 



for the spread of random walk, i.e., the mean square dis- 
placement ((x — x^Y) (X AT. The transient zero-energy 
manifold is given by if = (1/2)pq = 0. There is no 
stationary zero-energy manifold and no saddle point in 
accordance with the fact that random walk does not at- 
tain a stationary state. In Fig. El we show the simple 
random walk phase space. 



2. Overdamped oscillator 

The overdamped oscillator in ID is described by the 
Langevin equation 



dx 



-jx + r^it), (ritUO)) ^ AS{t), (A7) 



corresponding to the drift F ~ 2jx. We obtain the 
Hamiltonian 

H = - 7px = ^p{p - 2'yx), (A8) 



and ensuing equations of motion 
dx dp 





p 

» 








> 

E = 








< 




< 



(A9) 



E>0 



E>0 



FIG. 18: The {x,p) phase space for the random walk case. 
The transient zero-energy manifold is given by p = 0. 



with solution from a;' to x in time T, 



x{t) = ^^inh7t + x-sinh7(T-t) ^ ^^^^^ 



p{t) = 7 



sinh 7T 



sinh jt 



(All) 



From Eqs. (|A9p we have d^x/dt^ — 7^0; and the over- 
damped noise-driven oscillator thus corresponds to the 
motion of a particle in an inverted harmonic potential. 

The zero-energy manifold H = is composed of the 
transient submanifold p = and the stationary subman- 
ifold p ~ 2"fx intersecting at the saddle point {x,p) — 
(0,0). From the action S = (1/2) dt p{tY we obtain 
the familiar distribution 82] 



P{x, T) oc exp 



7 {x 



A 1 - e-27T 
For T — !■ 00 we infer the stationary distribution 



(A12) 



Po{x) (X exp 



7X 



(A13) 



In Fig. El we have shown the phase space for the over- 
damped oscillator. We note that in both the random 
walk case and the overdamped oscillator case, the weak 
noise method gives the exact and familiar results; quite 
similar to the corresponding calculations in the WKB 
approximation in quantum mechanics. 
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FIG. 19: The (x,p) phase space for the overdamped oscillator 
case. The transient zero-energy manifold is given by p = 
and the stationary zero-energy manifold by p = 27a;. The 
manifolds intersect at the saddle point (SP) (a;,p) — (0,0). 
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